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Abstract. 

We summarize the results of a recent study of heavy quark production 
and attenuation in cold nuclear matter. In p+p collisions, we investigate 
the relative contribution of partonic sub-processes to D meson production and 
D meson-triggered inclusive di-hadrons to lowest order in perturbative QCD. 
While gluon fusion dominates the creation of large angle DD pairs, charm on 
light parton scattering determines the yield of single inclusive D mesons. The 
distinctly different non-perturbative fragmentation of c quarks into D mesons 
versus the fragmentation of quarks and gluons into light hadrons results in a 
strong transverse momentum dependence of anticharm content of the away-side 
charm-triggered jet. In p+A reactions, we calculate and resum the coherent 
nuclear-enhanced power corrections from the final-state partonic scattering in 
the medium. We find that single and double inclusive open charm production can 
be suppressed as much as the yield of neutral pions from dynamical high-twist 
shadowing. Effects of energy loss in p+A collisions are also investigated in the 
incoherent Bertsch-Gunion limit and may lead to significantly weaker transverse 
momentum dependence of the nuclear attenuation. 
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1. Introduction 

A useful probe of the dense nuclear matter created in collisions of heavy nuclei 
at the relativistic heavy ion collider (RHIC) is one that is sensitive to dynamical 
scales and can be both cleanly measured experimentally and reliably calculated 
theoretically pQ. Because of color confinement, only hard probes, i.e. those with 
large momentum transfers, can be reliably calculated in the perturbation theory of 
Quantum Chromodynamics (QCD). Experimental measurements of inclusive particle 
suppression are now able to test jet quenching theory out to transverse momenta as 
large as pr ~ 20 GeV ■ On the other hand, typical dynamical scales of the nuclear 
matter produced in relativistic heavy ion collisions are on the order of hundreds of 
MeV, which is both much smaller than the scale of a hard probe and non-perturbative. 
Therefore, an ideal probe should be not only "hard" but also sensitive to this "soft" 
physics. Open charm production has a potential to satisfy these criteria because of the 
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two distinctive scales of the open charm meson: the charm quark mass, m c ~ 1.5 GeV, 
a relatively hard scale, and the binding energy, ~ Mjj — m c ~ hundreds MeV, which 
may be relevant to the fragmentation and dissociation of D mesons. 

With this in mind, we present a baseline study of heavy quark production and the 
first results on heavy-quark-triggered large-angle correlations in elementary nucleon- 
nucleon reactions p^. Next, we investigate the many-body QCD effects in single- 
and double-inclusive open charm production. Because of the quantitative importance 
of multiple scattering for inferring properties of the nuclear matter, the models of 
in-medium interactions should be subjected to experimental verification for heavy- 
flavor production. One may begin to gain confidence by applying the theoretical 
description of such effects to heavy ion collisions in simpler situations, such as proton- 
nucleus (p+A) scattering. In this case, plasma properties are not an issue and the 
main medium effect is, therefore, due to the multiple interactions in cold nuclei. The 
medium effects that we plan to investigate are also present in nucleus-nucleus (A+A) 
reactions during the interaction time Tint. = %Ra/i <C T eq . <C tqgp and cannot be 
neglected. Here Ra is the nuclear radius, 7 is the Lorentz gamma factor of the nucleus, 
T eg . is the equilibration time and tqgp is the lifetime of the plasma. 

We note that full technical details for the results presented in these proceedings 
are given in our complete manuscript pQ. 



2. D meson production and correlations in lowest order PQCD 

A formal expansion for the differential cross section of heavy quark production can be 
written as follows: 



da 



dyd 2 pT 
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where G(m,pr) —* 1 when m/pT — > 0, and power suppressed terms ~ (m/pr) n arc 
not shown. In Eq. JJJ A,B,a,i,bi are coefficient functions and logarithms arise due 
to the new mass scale, m. These are known to next-to-leading logarithm (NLL), but 
numerically, their contribution to the cross section is small up to pj> ~ 50 GeV 
Existing calculations of charm and bottom production [3] treat quarks as heavy, 
assuming <f> c ,b/N( x i A*/) = 0- We here evaluate the differential cross sections for 
open charm production and open-charm-triggered di-hadron correlations in the Born 
approximation. We include the charm contribution from the nucleon wave-function 
explicitly |T|, since this approach leads to a faster convergence of the perturbation 
series H] and correspondingly smaller next-to-leading (NLO) order K-factors. 

There is a direct relation between the absolute normalization of the cross section 
and the hardness of the non-perturbative D meson fragmentation. We compare in the 
left panel of Fig. ^ the calculated D° and D + cross sections to the Tevatron Run II 
y/s = 1.96 TeV data We find that three different combinations (Kjvlo — 1.5, 
r = 0.1), (Knlo = 2, r = 0.2) and (K NLO = 3, r = 0.4) yield little difference in 
the single inclusive charm meson spectrum. We note that for the same choice of the 
fragmentation parameter r — 0.1, a LO calculation with standard 4> c /n(%> Hf) 7^ 
gives open charm cross sections similar to the ones from a NLO calculation that treats 
flavor as "heavy" . 
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Figure 1. Left panel from 0: cross section for D° and D+ charm meson 
production to LO in perturbative QCD at t/s = 1.96 TeV. Data is from CDF 
Right panel from 1 : contribution of light hadrons and anticharm mesons to the 
fractional production cross section for pxi = 5 GeV D° + D + meson-triggered 
away-side correlations for yi = y2 = 0. The insert illustrates the difference in the 
fragmentation of light quarks and gluons into hadrons versus that for the c quark 
into charm mesons. 



Besides the inclusive spectra, we show our first results on charm meson-triggered 
away-side di-hadrons. The expectation for non-trivial pt 2 dependence of such large- 
angle correlations is based on the very different behavior with respect to z = 
phadron^pparton Q £ fragmentation functions of partons into light hadrons compared 
to those of heavy quarks into charm and beauty meson P^. The insert in the right 
panel of Fig.^shows that while light hadrons favor soft decays of their parent partons, 
heavy quark fragmentation is very hard. For z ~ 0.6 — 0.9 there is more than an order 
of magnitude enhancement in the decay probabilities Ddo/ c (z) + Dd+/ c {z) relative 
to (l/N)J2i = i jv(-^w±/i + Dfr±/i + D p rp\u), where i runs over the light and charm 
quarks, antiquarks and the gluon. 

Triggering on a D meson fixes the momentum of the charm quark much more 
reliably than does triggering on a light hadron |S] in studying away-side correlations. 
The non-perturbative fragmentation, therefore, controls the yields of di-hadrons versus 
the associated momentum px 2 and the abundances of the different particle species. 
The right panel of Fig. ^ shows our prediction for the hadronic composition of the 
D° + D + -triggered away-side jet. At transverse momenta significantly smaller than 
the trigger transverse momentum, pt 2 <C pr t , the away-side jet is dominated by pions, 
kaons and protons. At transverse momenta ~ Pt 2 -, the away-side jet is expected to 
be dominated almost completely by D° and D~ mesons. Since there is little sensitivity 
to the choice of r, our prediction is robust and can be used as a test of the production 
mechanism of heavy quarks. 

3. Resummed QCD power corrections to open charm production 

In the case of hadronic collisions there is indeed similarity with the DIS dynamics in 
the final-state rescattering of the struck, small X},, parton from the nucleus, as shown 
in the left panel of Fig. [21 Such nuclear-enhanced power corrections are equivalent 
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to dynamical mass for the interacting quark, similar to the generation of dynamical 
electron mass in a strong magnetic field, 
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They lead to suppression of the single- and double-inclusive hadron production cross 
sections in deep inelastic scattering (DIS) [7] and forward rapidity p+A reactions 
as long as the minimum coherence criterion x < 0.1 is satisfied. 

The mechanism behind the reduction of the deep inelastic inclusive scattering 
cross section and the differential hadron production cross sections is illustrated 
schematically in the right panel of Fig. [21 PDFs are falling functions of the momentum 
fraction x. The evaluation of the parton distribution function at an effectively larger 
value of x leads to a reduction of the cross section. From our theory we can easily 
derive a hierarchy of dynamical high-twist shadowing, 

-Act(.t, Q 2 

for cross sections with multiple scattering dominated by valence quarks, sea quarks or 
gluons, respectively. 

In £Q, we completed the summation of all- twist, nuclear-enhanced power 
corrections for massive-quark final states. In hadronic collisions, we identify the 
function that contains the small Xb dependence F a b^ c d(xb) = 4>b/N(xb)\M a b^ c d\ 2 /xb- 
Implementing the normalization to no nuclear effect on a single nucleon as follows, 
A 1 / 3, — * A 1 / 3, — 1, and having allowed for both massless and massive partons (quarks 
or gluons) we find pQ 

^(A 1 / 3 — 1) 



S(x,Q 2 ) = 



S v (x,Q 2 ) < S s (x,Q 2 ) < S g (x,Q 2 ) 



(3) 
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In Eq. Cd = 1 (9/4) for quarks (gluons), respectively. When m d — ► we recover 
the result of 6 . The mass here reduces the effect of power corrections at — i ~ m 2 ,. 

While we have focused in detail on the power corrections in the i-channel, it 
should be noted that the coherent initial- and final-state interactions of the incoming 
parton from the proton have also been calculated pQ. The results, that we present 
here, are analogous to Eq. J3J, 



Fab^cd(xb) => F a b^cd Xb 



F a b->cd(xb) => F a b^cd X b 
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By direct comparison of the nuclear-A 1 / 3 -enhanced power corrections in Eqs. I@J), (0 
and © it is easy to see that in the forward rapidity region \t\ <C \s\, only DIS-likc 
dynamical shadowing corrections are important. 

In processes with no strong final-state interactions, such as the Drell-Yan, 
coherent multiple initial-state scattering can give only enhancement of the cross 
section, as shown in Eq. JSJ, because of the rescaling of Xb to smaller versus larger 
values. Thus, initial-state multiple soft interactions without energy loss always lead 
to a Cronin effect, independently of whether they are treated as coherent or not. The 
derivations, summarized here, are critical to elucidating the dynamical origin of nuclear 
effects in high energy hadronic reactions. We have shown, through explicit calculation, 




Figure 2. Left panel from 0: multiple coherent scattering of the outgoing 
partons in proton-nucleus reactions in the t-channel. We have denoted by -x-> and 
— |— » the long distance and contact propagators, respectively. We have indicated 
the lightcone positions, e.g. yf . Arcs show the momentum routing. Right panel 
from 0: effects of rescaling of the momentum fraction x illustrated on the example 
of the lowest order CTEQ6 parton distribution functions for valence and sea 
u-quarks. Gluons rescatter in the final-state to lowest order only in hadronic 
collisions and are included for completeness. 



that such effects are process dependent and may change their sign. Therefore, these 
are not factorizable as a part of the parton distribution functions and fragmentation 
functions. 



4. Energy loss in cold nuclear matter 

In the presence of a nucleus, final-state rescattering of the struck small- x parton 
with its remnants exhausts the similarity between hadronic collisions, such as p+A 
and A+A, and DIS [Sj. In I + A (DIS), the multiple interactions of the incoming 
leptons are suppressed by powers of a em /ct s relative to the struck parton scattering. 
In contrast, in p+A and A+A the initial- and final-state scattering of the incoming 
quarks and gluons are equally strong. Nuclear modification, in particular jet energy 
loss associated with the suppression of particle production, cannot be neglected. 

To illustrate the importance of energy loss, we study the nuclear modification 
of hadron production over a large range of center of mass energies and momentum 
fractions x^. We implement the presently well known incoherent limit of energy loss 
for asymptotic t = — oo on-shell partons first derived in [Sj. The double differential 
gluon intensity spectrum per scattering and the fractional energy loss e read: 

W 2 k\ 2 k 2 (k-q) 2 ' £ A ' {) 

where k is the transverse momentum of the radiative gluon, q is the transverse 
momentum transfer from the medium and u> is the gluon energy. In Eq. |7| l L/A 
is the number of interactions and k = 0.0175 for minimum bias reactions implies that 
an average parton loses ~ 10% of its energy in a large nucleus such as Au or Pb. 

We first examine the results obtained by the CERN NA35 fixed target experiment 
with y cm = 3. It measured hadron production in d+Au reactions at ^Jsnn = 
19.4 GeV JOI- We take R v (pt), the ratio of hadron spectra in different rapidity 
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Figure 3. Left panel from 0: nuclear modification at three different rapidities 
y — Vcm in y/sjvjv = 19-4 GeV d+Au collisions. Data is from NA35 1101 . Right 
panel from 0: calculations of dynamical shadowing, with and without cold 
nuclear matter energy loss, are compared to the measured fill 7r° suppression 
at y = 4 for d+Au collisions at y'sjvlv = 200 GeV. 



bins, and Raa(pt), 

R AB (p T ) = da AB(y,PT)/dyd 2 PT Raa (pt ) = dN AB(y,PT)/dyd 2 PT , g , 
V da\ B (y ha se,PT)/dyd 2 pT ' T A B(b)da% N (y,p T )/dyd 2 p T ' 

as measures of nuclear matter effects. The left panel of Fig. [3] shows this forward 
rapidity suppression of negative hadrons relative to the baseline y — y cm rs —1.5 
production cross section. Dynamical shadowing calculations 6] give < 5% effect at 
this energy in this rapidity range. Conversely, the implementation of energy loss in cold 
nuclear matter leads to much larger suppression at forward rapidity and significant 
improvement in the theoretical description of the data in Fig. |3 

In the other extreme, at y/s — 200 GeV d+Au collisions at RHIC, the STAR 
collaboration observed a factor of 3 suppression of 7r° production at rapidity y = 4 . 
The right panel of Fig. [3] shows that dynamical high-twist shadowing, constrained 
by the DIS data down to values of Bjorken xb ~ 10~ 4 , underpredicts the nuclear 
attenuation by a factor of 2. Incorporating the same energy loss that we identified at 
much smaller center of mass energies, we observe very good agreement between theory 
and the data [Q. 



5. Numerical results 



To illustrate the similarities and differences between massless and massive final-state 
partons, we carry out a comparative study of the effect of power corrections on single 
and double inclusive 7r° and D meson production. The left panel of Fig. 0] shows 
the suppression of the low and moderate p? 1 neutral pion cross section from high- 
twist shadowing relative to the binary scaled p+p result. The nuclear modification 
factor is shown for two different forward rapidities, y\ = 1.25,2.5, and three different 
centralities, b — 3 fm, 6 m in.bias = 5.6 fm, b — 6.8 fm, in ^/s = 200 GeV d+Au collisions 
at RHIC. At transverse momenta pr t = 1.5 GeV the suppression can be substantial, 
but disappears toward higher pr l due to the power law nature of the effect. Coherent 
power corrections cannot fully account for the nuclear suppression measured by the 
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Figure 4. Left panel from 0: suppression of the single inclusive hadron 
production rates in d+Au reactions versus pT t for rapidities y± = 1.25 (smaller 
attenuation) and y± = 2.5 (larger attenuation). Theory does not include energy 
loss. Right panel from J]: the same calculation with cold nuclear matter energy 
loss of the Brtsch-Gunion type |S] included. Data is from PHENIX Q^J IT3I . 
Bottom panels show impact parameter dependence of the calculated nuclear 
modification for central, 6 = 3 fm, minimum bias, b m i n .bias = 5.6 fm and 
peripheral, b = 6.8 fm, collisions. 



PHENIX experiment J2| an d the discrepancy becomes larger at higher transverse 
momenta. The calculated suppression of D° + D + mesons (and equivalently, D° + D~ 
mesons), in deuteron-gold collisions at RHIC is also shown in the left panel of Fig.0J It 
should be noted that the nuclear modification is very similar to or slightly larger than 
that for 7T°. The reason for this similarity is that in both cases the dominant channel 
of hadron production is via quarks scattering on gluons, which have the same (large) 
color singlet coupling to the medium, Cd = 9/4 in Eq. Q). In addition, the typical 
momentum fraction, z\, which enters in the determination of the hard scale, t oc X/zi, 
is slightly larger for the D mesons. Preliminary PHENIX data on the modification of 
ji~ coining from the decay of heavy flavor at y = 1.6 in cold nuclear matter is provided 
as a reference ^21- Although the error bars are large, D mesons seem to be suppressed 
as much as light hadrons. 

Having investigated the largest dynamic range of measurements in proton(deuteron)- 
nucleus reactions accessible to perturbative QCD calculations in the previous section, 
we return to D meson production and correlations at RHIC. We use the same frac- 
tional energy loss e = AE/E as in Fig. 01 Single inclusive tt° and D +D + suppression 
at rapidities y\ — 1.25, 2.5 and minimum bias, central and peripheral d+Au collisions 
at RHIC are shown in the right panel of Fig. In this case, very good agreement 
between the QCD theory incorporating cold nuclear matter effects and the PHENIX 
measurement of muons coming from the decay of light hadrons ^21 is achieved. We 
find that the magnitude of the D meson suppression is similar to that of pions. It is 
also comparable with the first forward rapidity results on heavy quark nuclear modifi- 
cation in d+Au reactions at RHIC ^2] ■ Nuclear modification of inclusive two particle 
production is seen to follow closely that of single inclusive hadrons pp. 

There are similarities between the calculations that only include resummed 
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nuclear-enhanced power corrections, left panel of Fig. 0] and the calculations that do 
not neglect energy loss in cold nuclear matter, right panel of Fig. 21 Both effects are 
generated through multiple parton scattering in the medium and lead to suppression 
of the rate of hard scattering P^. In both cases single inclusive particle production 
and large angle di-hadron correlations are similarly attenuated. Like all nuclear many 
body effects, these increase with the centrality of the collision. The difference in the 
resulting nuclear modification is that high-twist shadowing arises from the coherent 
final-state scattering of the struck small-Xfc parton of the nucleus and disappears 
as a function of the transverse momentum |S] . The energy loss considered here arises 
from the initial-state inelastic scattering of the incoming large- x a parton from the 
proton(deuteron) |5] and leads to a suppression which is much more px independent, 
similar to final state energy loss applications 0- 



6. Conclusions 



In these proceedings, we summarized the results of the first perturbative QCD 
calculation of heavy meson triggered large-angle di-hadron yields pQ and showed how 
such measurements can provide information on the D meson production mechanism 
and the non-perturbative fragmentation of heavy quarks. We identified multiple 
scattering effects that drive the cold nuclear matter attenuation of D meson and 
light hadron production at forward rapidities, namely, high twist shadowing |SJ [7] 
and parton energy loss (Sj 03 • Our work demonstrates P how such effects 
can be independently constrained by DIS data (shadowing) and low energy p+A 
measurements jl(J| where coherent power corrections are not important (energy loss). 
In the framework of the established QCD collinear factorization, our approach of 
incorporating dynamically generated and process-dependent nuclear effects is expected 
to give a more reliable description of particle production ^2 El H3| m P+A collisions. 
This work also provides the baseline for precision QGP tomography using heavy quarks 
as probes of the plasma away from midrapidity pQ where new experimental capabilities 
at RHIC and the LHC are expected to become available. 
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